DELTA LAUNCH VEHICLE ISOGRID STRUCTURE NASTRAN ANALYSIS 


By Daniel J. Knighton 
NASA Goddard Space Flight Center 


SUMMARY 


The purpose of this paper is to present the varied structural analysis 
applications of NASTRAN on the new isogrid 1 structure of the Delta launch 
vehicle. Isogrid .is a method employed to stiffen the booster cylindrical shell 
surface by integrally machining ribs and skin from one piece of metal. The ribs 
are arranged in a repetitive equilateral triangular pattern; this results in a 
structural surface whose stiffness is orthogonally isotropic. For that struc- 
ture both static and general instability buckling analysis will be described 
and the theoretical results will be compared to test data. One of the param- 
eters that has been included in this analysis is the effect of pressure load- 
ings on the buckling allowable load. In addition, a description of a supple- 
mentary GSFC -originated computer program, without which this analysis could not 
have been performed, will be presented. 


INTRODUCTION 


A brief description of the Delta launch vehicle and its performance 
advancements will be given with specific attention focused on the IBM 360/95 
computer analysis of the recently innovated Delta external isogrid shell struc- 
ture. The NASA STructural ANalysis (NASTRAN) program has been successfully 
employed to (l) verify the isotropic property of isogrid, (2) show the variation 
in deflections and detail stress levels depending on loading versus the orienta- 
tion of isogrid, (3) perform as an engineering liaison tool, and (4) determine 
the general instability buckling allowables for isogrid tank structures. 


-^-Isogrid structure concept was developed for Delta by McDonnell Douglas 
in 1970. 
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side of isogrid triangle, cm (in.) 

area, cm 2 (in.) 

locations of CBAR 

NASTRAN bar element 

NASTRAN triangular plate element 

rib height, cm (in.) 

diameter of cylinder, cm (in.) 

Young's modulus of elasticity, kN/m 2 (psi) 
isogrid modulus, kN/m 2 (psi) 
force, N ( lb) 

rib-allowable compression stress, kN/m 2 (psi) 
skin-allowable compression stress, kN/m 2 (psi) 
force, N ( lb) 

modulus of rigidity, kN/m 2 (psi) 
height of isogrid triangle, cm (in.) 
moment of inertia, cm 4 (in. 4 ) 
polar moment of inertia, cm 4 (in. 4 ) 
form factor of cross section 
elastic stiffness 
differential stiffness 
length, cm (in.) 
moment, N-m (in. -lb) 
bending moment, N-m (in. -lb) 

isogrid general instability allowable load, kN/m (ib/in.) 
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rib-allowable compression or crippling load, klf/m (ib/in.) 

skin- allowable compression or buckling load, kN/m (lb/in.) 

equivalent compression load, kN/m (lb/in.) 

pressure, kW/m 2 (psi) 

applied axial load, N (lb) 

radius, cm ( in. ) 

equivalent skin thickness, cm (in.) 
rib thickness, cm (in.) 
skin thickness, cm (in.) 

weight equivalent thickness (accounting for weld land and/or edge 
panel design) 

torque, N-m (in. -lb) 

linear displacement of CBAR grid point 

eigenvector 

shear, N ( lb) 

infinitesimal element volume 
extensLonal strain, cm/cm (in. /in.) 
shear strain, cm/cm (in. /in.) 
rotation of CBAR grid point 
eigenvalue 
Poisson’s ratio 

axial stress of isogrid in x, y directions, kN/m 2 (psi) 
shear stress, kN/m 2 (psi) 
rotation of dV 
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DELTA ISOGRID STRUCTURE 


The Delta launch vehicle project began in 1959 shortly after the forma- 
tion of MSA. It started with existing stages that were comprised of the Thor 
first stage and Vanguard second and third stages. Since that time it has been 
project policy to upgrade the Delta booster for increased performance and 
reliability at a minimum cost growth as shown on figure 1. The current launch 
vehicle and improved configurations are shown on figure 2. 

One of the many upgraded design features of these new Delta vehicles is 
the incorporation of a cylindrically shaped isogrid structure. This design 
forms the external skin structure of the booster. The new tank skin is 
initially machined from flat 1.27-cm (l/2-in.) thick 14ST 6 free -machining 
aluminum alloy plate, brake-formed into curved shapes, and finally welded into 
2.44 -m (8-ft) diameter tank shells. Typical integrally machined structures, 
waffle and isogrid, are shown in figure 3* Isogrid now replaces the less 
efficient rectangular waffle tank and skin, stringer, and frame interstage and 
fairing structures. The improved Delta isogrid tank construction is about 50 
percent more efficient (weightwise), than the waffle and equivalent in efficiency 
to the built-up assemblies with regards to compressive local and general buck- 
ling capabilities. It is also interesting to note that isogrid is orthogonally 
isotropic in its overall strain characteristics. 

The first vehicles to employ the isogrid structure are now planned to 
carry the IMP H and TELESAT A spacecraft into orbit by the end of this year. 
Isogrid hardware is not only being fabricated for Delta vehicles but it is also 
being developed for the Shuttle (ref. l) and Skylab programs under the cogni- 
zance of MSFC and is being proposed for more efficient Delta/Agena replacement 
second stages (fig. 4) in the Versatile Upper Stage (VUS) study (ref. 2) to 
GSFC requirements. Preliminary evaluations performed on Delta and the VUS in- 
dicate that isogrid structures are approximately 50 percent the cost of skin, 
stringer, and frame structures. This was the reason why the interstage and 
fairing were designed using isogrid for new model Deltas. 


ANALYSIS OF ISOGRID BY MSTRAN 


The new structure, a relatively fine rib-node mesh, was modeled using 
NASTRAN program elements and procedures. Both static and buckling analysis 
were performed to give the Delta project a good understanding of the new struc- 
ture being designed for Delta. Loadings applied to the isogrid model initially 
were only compression type but later were expanded to include pressure effects; 
i.e., hoop loadings. The compression loading included axial compression, equiv- 
alent axial compression due to bending, and offsetting axial tension forces 
pR/2 due to internal pressure,. The buekling analysis also included, for com- 
parison only, the effects of an external pressure loading. 
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PROOF OF ISOTROPY OF ISOGRID 


For an isogrid orientation having 
the &' and \f5&' periodic grid distances 
shown to the right, it has been proven 
(ref. 3 ) that the material has isotropic 
elastic properties and that it obeys 
Hooke ' s law . 
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where the isogrid elastic modulus E is related to the actual material modulus 
E by E = t r E/h. 


Subsequent to the previous analysis, the rib__and plate elements were 
combined to produce an equivalent skin thickness t such that the booster tank 
could be treated like an equivalent monocoque cylindrical structure. The 
equivalent monocoque skin thickness is 


t = 



t-£»d 

IT 
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1.270 cm 
(0.500 in.) 


The t relation is important in solving for skin and rib compression allowables 
(ref. 4), which are 

Skin buckling : 

N sb = 10.875Et(^) 

Esb 
^sb = T~ 


125 



Rib crippling: 


N rc = 0.422Et 

Nrc 
frc = -f 

Isogrid, being isotropic, does possess an idiosyncrasy, i.e., its detail 
stresses must vary internally within the isogrid depending on the direction in 
which loadings are applied. The reason for this is that the isosceles tri- 
angular arrangement is oriented slightly differently in the x- and y-axis direc- 
tions. To learn more about this novel structure, a small finite element model 
was fashioned using NASTRAN bar (CBAR) and triangular plate (CTRIA2) elements, 
figure 5 . Specifically, this model was prepared to check the overall isotropic 
property of the integral rib and skin combination and to review the internal 
rib and skin stresses that result when the same magnitude of loading is 
placed in either the x or y directions. 

The compression loading was selected from a typical Delta loading profile 
and is equated as follows: 

N x = ~ + = 308.2 kN/m (1760 lb/in.) 

TrD 1 tu 

N x was applied first in the negative x and then in the negative y direc- 
tions. The results, including deflections and stress levels, are shown in 
figures 6 and 7» Note that the overall deflections, 1.346 mm ( 0.0530 in.), are 
the same for either x- or y-axis loadings, which indicates the overall iso- 
tropic characteristic of isogrid. Further, the rib and skin compressive stress 
levels in the center bay were lower for the x-axis loading direction (loads 
perpendicular to the main ribs) by 31 and 1 percent, respectively. Thus, to 
keep the compressive stresses at a minimum in the isogrid propellant tanks, the 
main ribs were located accordingly (on Delta stations, fig. 5). 

The initial model, figure 5 , was 46.94 cm ( 18.48 in.) square and required 
additional plate and bar elements. Two computer runs were made using heavy 
load support elements (50 times stiffer than isogrid CBAR members) and multiple 
point constraints at the load grid points and the same deflections and stress 
levels were obtained. The grid-point restraints are as noted on figure 6. 

A third computer run was accomplished using a simpler model as shown in 
figure 8 . Note that the model is 46.94 by 50.80 cm (l8.48 by 20.00 in.) and 
contains fewer plate elements. The deflection and stresses are as shown in 
figures 9 and 10. The deflection and stresses are approximately the same as 
produced by the previous more complicated model. This simpler model was also 
utilized for the combined loading analysis, which is described in the next 
section. 
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COMBINED LOADING OF ISOGRID 


The previous NASTRAN models were loaded in compression only. Next, the 
preceding model, having a single fixed point (grid point 28), was loaded with 
axial compression (308.2 kN/m (l?60 lb/in.)) and hoop tension (4l4 kN/m 2 (60 
psi)) forces. For the levels of loading applied, the results are as shown in 
figures 11 and 12. The rib compressive stress -level was 51 percent lower and 
the skin buckling level was 3 percent higher for the x-axis loading conditions 
where the compression loadings were perpendicular to the main ribs (fig. ll) . 

A close look at the actual rib and skin stresses and a direct comparison with 
their respective allowable stresses is given in table 1. Note that the skin 
actual stress was closer to its allowable stress. For the combined loads case 
analyzed, the skin stress will be a minimum if the compression load is in 
line with the isogrid main rib. Because past tests (ref. 5) indicate that 
local skin buckling will reduce general instability allowables (by approximately 
20 percent), it is important to keep the skin actual stresses as low as pos- 
sible. 


TABLE l.--Rib and Skin Compression Stresses 


Element 

Actual 
stress, 
kN/m 2 (psi) 

Allowable 
stress, 
kN/m 2 (psi) 

/Allowable stress -A -inn 

\ Actual stress / ’ 

percent 

Rib 

-252 067 
(-36 558) 

-424 042 
(14ST651 yield 
corrected for 
plasticity; 
-61,500) 

+68 

Skin 

-173 637 
(-25 183) 

-202 024 
( -29 300) 

+l6 


LIAISON COMPUTER ANALYSIS 


During the early stages of fabrication, there have been occasions when 
isogrid skins have been machined improperly. In the cases where skin pockets 
and rib thicknesses were undersized at random locations, the small NASTRAN 
model previously described was used effectively. Having the actual stress 
levels throughout the model, it was quite easy to determine rib and skin margin 
of safety and whether reinforcements were required. 


127 









ISOGRID TANK STATIC AMD BUCKLING ANALYSIS 
(MOVE AND NASTRAN PROGRAMS) 


This analysis was performed to determine the general instability theo- 
retical allowables for axial compression and various tank hydrostatic loadings; 
i*e., negative or external, zero, or positive or internal pressure loadings. 

The fact that tank pressure loading can be handled quite easily by the 
NASTRAN program made this a more interesting and comprehensive study. The 
results of this analysis compared with test data and other methods of analysis 
(ref. 6) will be summarized later. Before discussing the model and the results 
a brief description of the theory will be presented. 


BUCKLING THEORY 


NASTRAN buckling analysis (ref. 7 ) is based on the elastic and differen- 
tial stiffnesses of the structure analyzed. The elastic properties of a struc- 
ture are generally dependent on shear AG/K, torsion JG, bending El, and 
axial AE stiffness characteristics. The differential stiffness is based on 
the static loading, displacement, and geometry of the structure. An example of 
a bar element differential stiffness matrix is shown in figure 13 (ref. 7)* 

The approach presented in reference 7 is essentially based on using Lagrange's 
equations of motion on a structural system with a finite number of degrees of 
freedom. 

The steps (ref. 7) for including differential stiffness in a structural 
problem are as follows: 

(1) Solve the linear static response problem for the structure in the 
absence of differential stiffness and compute the internal forces 
in each element. 

(2) Using the results of (l), calculate the differential stiffness 
matrix for the individual elements and apply the standard reduction- 
procedure (constraints and partitioning) to form the differential 
stiffness matrix [Kaa d l in final form. 

( 3 ) For the buckling problem, find the characteristic eigenvalue and 
eigenvectors for 

l*aa + AK aa d ] (u) = 0 

For the isogrid cylinder buckling problem, the inverse power method was used to 
determine the eigenvalues. Because the above-defined set of buckling equations 
are homogeneous, the assignment of one arbitrary eigenvector (within NASTRAN) 
must be made to start the analysis. The solution of the first trial eigenvalue 
permits the determination for convergence; the redetermination of the eigen- 
value, if necessary; and the subsequent final solution of all eigenvector 
ratios {u} or of the mode shape. 
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It is the eigenvalue that is the factor by which the applied compression 
static loading is multiplied to produce the theoretical threshold of buckling. 
Once obtained, the theoretical buckling load is multiplied by a 65 -percent 
reduction factor. That factor has evolved from testing full-scale isogrid 
cylinders and it accounts for manufacturing tolerances and distortions that 
are present in hardware. 


BUCKLING MODEL DESCRIPTION 


During June 1971? three isogrid tank segments (244-cm diameter by 244-cm 
length (8 -ft diameter by 8 -ft length)) were fabricated from free-machining 
2014 t 6 aluminum; subsequently, compression tested with zero internal pressure. 

The mathematical model was determined, based on those compression test tank 
sizes, in order to correlate computer results with test values. 

The NASTRAN model, being an exact replica of the tested tank segments, 
was a rather fine and detailed model. This model size is believed to be the 
largest structural model analyzed at NASA Goddard Space Flight Center. The 
model was composed of NASTRAN triangular plate (CTRIA2) and bar (CBAR) elements. 

A complementary computer program, MOVE (ref. 8), was utilized to simplify the 
task. A model such as ours, having 1850 grid points and 7986 elements, would 
have been extremely difficult, if not impossible within the time available, to 
construct had it not been for the help of the MOVE program. 

The MOVE program was designed to generate repeat bulk data for structures 
having a number of identical segments. All one need do is to generate one 
basic NASTRAN bulk data segment and MOVE does the rest. The strip model 
selected for our problem was a longitudinal segment shown in figure 14. By 
moving that segment 66 times (5-45° per move circumferentially), the complete 
cylinder was generated as shown in figure 15 . It should be noted here that all 
of the undeformed and deformed figures in this report were made by the Stromberg- 
Carlson 4o60 plotter, which is one of several systems that may be specified in 
NASTRAN programing. 

When the MOVE program was first programed, it was designed for a relatively 
small structural problem. Thus, the initial trial of that program on a large 
structural model was not successful. The reason was that the MOVE program 
could not handle CBAR continuation cards and a very large number of repeat 
cycles. Further, the MOVE program was written to place its bulk data output 
on punched cards. This feature could have been quite cumbersome to handle for 
our large buckling problem ( approximately 8 data card file boxes, or 8000 
cards). Including the output data on tape was a far superior way of handling 
an immense quantity of data. That program, written in PL-1 language, was 
subsequently revised to rectify the above problem areas and to store its out- 
put data on magnetic tape. 

Having accomplished the MOVE programing successfully, the next step was 
to input the data using NASTRAN programing methods. Here too, initial attempts 
were not good because of an assortment of program problems involving the 
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application of the MOVE data tape and the determination of a practical size 
storage requirement for our isogrid structural model. The solutions were 
readily available by properly arranging the executive control section of 
RASTRAN and adjusting the storage requirement from 500 000 to 1 000 000 places. 
Incidentally, the isogrid tank buckling computer analysis required all but 5000 
places for some of the RASTRAN routines --that was measuring storage quite close. 
Yet still another problem occurred: The computer program timed out at 120 min 

of requested central processing unit (CPU) time. Rather than increase the run 
time, the model was reduced, but reduced in such a way that the original size 
tank configuration could be maintained. This was achieved by reducing the 
RASTRAN model to one-quarter the original size, limiting the buckling problem 
by analyzing only the first four symmetric-symmetric modes as shown in figure 
1 6, and establishing the constraints or degrees of freedom accordingly. 

To accomplish the reduction in model size to one-quarter the original, 
the grid points were modified to include 1 through 13 (fig. 14) only and the 
MOVE rotations ( 5 .45°/rotation) were reduced from 66 to 33* The MOVE program 
was rerun and the reduced data ( on tape) were submitted via RASTRAN for a plot 
of the one -quarter -scale model shown in figure 17. 

The buckling analysis was again attempted and constraint problems were en- 
countered. A more careful look at the symmetric-symmetric modes and the loaded 
end conditions led to grid-point degrees of freedom revisions that subsequently 
produced good buckling results as shown in table 2. Plots of the statically 
deformed and buckled structure are presented in figures 18 to 24. 

But how dependable was the general instability buckling analysis using 
NASTRAR? To answer this question, the RASTRAN output was compared directly to 
another analysis (ref. 6) and to isogrid tank compression test data (refs. 4 
and 9) as shown in figure 25. For that figure the theoretical buckling values 
were multiplied by a 65 -percent test correlation factor (ref. 4) and compared 
to 99 percent probability buckling test values (2.34a lower than the average of 
nine test data points). The lower than theoretically predicted buckling allow- 
ables thus achieved account for imperfections that are inherent in large struc- 
tures. Such imperfections are generally attributable to built-in residual form- 
ing stresses, slight amounts of tank out -of -roundness, and local rib and skin 
waviness . 

Figure 25 graphically documents instability buckling levels for the Delta 
isogrid and waffle designs as compared to monocoque designs. It also shows the 
allowable skin buckling running load for the isogrid design. The graph has 
been produced by placing the weight equivalent thickness (t^t. eq.) along the 
abscissa and the allowable compressive load (Ngi or Rsb) along the ordinate. 

By scanning figure 25, it becomes apparent that the general instability 
allowables obtained through the use of NASTRAN correlate very well with the test 
data point and the cylindrical shell analysis by ¥. Flugge (ref. 6) and McDon- 
nell Douglas H312 analysis. Because the symmetric-symmetric buckling results 
appeared to be quite good, further analyses in the asymmetric buckling modes 
were not performed. 
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Note also that the instability plots are based on compression loadings and 
no internal pressure. To verify the pressure effect, analysis was performed 
using NASTRAN for the conditions of compression loading and 4l4-kN/m 2 (60-psi) 
internal and external pressure. These data points are also shown in figure 25. 
It is interesting to observe that internal pressure (through pressure stiffen- 
ing of skin and rib) improves the zero-pressure buckling allowable by 7 percent, 
and the external pressure degrades the zero-pressure buckling level by -22 per- 
cent. The 7 -percent improvement correlates fairly well with a 4 -percent gain 
seen during testing of scaled-down lexan isogrid models (ref. 5)* 

Observe also that isogrid is local skin buckling critical from a thickness 
of 2.39 mm (0.093 in.) and lower. During the previously mentioned scaled-down 
lexan isogrid test program, local skin buckling levels with and without internal 
pressure were evaluated. The results are what would be expected- -the general 
instability buckling load level was improved by 20 percent by adding an internal 
pressure loading. The reason for that improvement was that internal pressure 
stiffened the ribs and skin pocket areas such that 100 percent of the theoreti- 
cal general instability level could be achieved. It is for this reason that a 
combined loads (compression, shear, and internal pressure) research and develop- 
ment task has been proposed. The combined loads would be applied to several 
'full-scale aluminum isogrid tank segments (2.44-m diameter by 2.44-m length 
(8-ft diameter by 8-ft length)). 


CONCLUSIONS 


The NASTRAN program techniques have been advantageously employed on the 
Delta launch vehicle isogrid structure in the following areas: 

(1) The combined loads isogrid model served as an excellent analytical 
tool in accomplishing liaison stress analysis. That isogrid is 
orthogonally isotropic was proved by NASTRAN. 

(2) NASTRAN buckling analysis produces good general instability buckling 
allowable load levels. It also has the capability of combining axial 
compression, bending, shear, and pressure loadings. Shear effects 
were omitted in this analysis in order to correlate test and other 
analyses. 

( 3 ) Buckling analysis of the isogrid cylinder determined that internal 
pressure loading, when combined with axial compression loading, 
provides a 7 -percent improvement in the general instability allow- 
able. Because both local skin pocket and instability buckling 
allowables are improved by internal pressure stiffening effects, a 
combined loads development task has been proposed. 
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SPACECRAFT WEIGHT IN SYNCHRONOUS TRANSFER ORBIT (lb) 



Figure 1. — Delta performance and cost history 
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•Typical of thisvehicle and more advanced models. 


Figure 2. --Configuration development plan 
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£ DELTA STATION 
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1.27 mm 
(0.050 in.) 

Hh 


1.270 cm 
(0.500 in.) 
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1 j!7 mm 
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1.19 mm 
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Figure 5. --Delta isogrid NASTRAN 
model . 


COMPRESSION 
308.2 kN/m 
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TRIANGULAR PLATE 72T2 
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Figure 6. --Deflection and stresses, 
compression load in -x direction. 
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CTR1A2 TRIANGULAR PLATE ELEMENTS 


GRID POINTS 



Figure 14 . - -NASTRAN segment for MOVE program. 



TYPICAL ISOGRID SECTION 


Figure 15 - --Isogrid tank structural 
model (1850 grid points; 7986 ele- 
ments) , 
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CIRCUMFERENTIAL MODES 



MODE FREQUENCY, CYCLES: 2 



LONGITUDINAL MODES 


2.44 m (8 ft) LONG 


MODE FREQUENCY, CYCLES: 'h 




Figure 16. --Symmetric-symmetric modes for isogrid 
tank 2.44 m (8 ft) in diameter and 2.44 m (8 ft) 
in length. 


(£ 2.44-m (8 ft) LONG CYLINDER 



Figure 17. --One quarter tank model undeformed (l.2-m 
(4 -ft) length, 1.2 -m (4 -ft) radius). Axial compres- 
sion forces include bending moment equivalent force, 
axial force, and offsetting pressure equivalent 
force. 
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Figure 18. --Static de- 
flections, compression 
load plus 4l4-kN/m 2 
(60-psi) external 
pressure. 



Figure 19. --First 
buckling mode, com- 
pression load plus 
4l4-kN/m 2 (60-psi) 
external pressure. 



Figure 20. — Second 
buckling mode, com- 
pression load plus 

4l4-kN/m 2 (60-psi) 

external pressure. 



Figure 21. --Static de- 
flections, compres- 
sion load plus 4l4- 
kN/m 2 (60-psi) in- 
ternal pressure. 
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Figure 22. --First buckling 
mode, compression load 
plus 4l4-kN/m 2 (60-psi) 
internal pressure. 



Figure 23- --Static deforma 
tion, compression load 
plus zero internal 
pressure. 



Figure 24. --First buck- 
ling mode, compression 
load plus zero internal 
pressure. 
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or N . h , Ib/in. 



NASTRAN GENERAL INSTABILITY (Gl) ALLOWABLE, 

D = DIAMETER = 244 cm (96 in.), t s “ 1 .27 mm 
(0.050 in.) TO 1.65 mm (0.065 in.), AND t, = 1.78 mm 
(0.070 in.) TO 2.03 mm (0.080 in.) 

□ 414-kN/m 2 (60-psi) INTERNAL PRESSURE 
O ZEROpsi INTERNAL PRESSURE 
A 414-kN/m 2 (60-psi) EXTERNAL PRESSURE 

REFERENCE 6 ANALYSIS OR MDAC H312 PROGRAM 

SKIN BUCKLING ALLOWABLE DERIVED FROM 

ISOGRID TEST 

f DELTA ISOGRID TEST CYLINDER, 99 PERCENT 

PROBABLE LOW VALUE, D » DIAMETER = 244 cm 
(96 in.), t s = 1.65 mm (0.065 in.), t, = 2.03 mm (O.OBO in.) 
WAFFLE AND MONOCOQUEGI ALLOWABLES 
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Figure 25 . --Buckling allowables, 14-ST6 aluminum cylinders 
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